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SUMMARY 


A stability analysis is made of a long flat rectangular plate sub- 
jected to a \iniform longitudinal compressive stress and supported along 
its longitudinal edges eind along one or more other longitudinal lines 
by elastic line supports. The elastic supports possess deflectional 
and rotational stiffness . The results of the analysis are presented 
in the form of charts in which the buckling-stress coefficient is plotted 
against the buckle length of the plate for a wide range of support stiff- 
nesses . The charts make possible the determination of the compressive 
buckling stress of plates supported by members whose sniffness may or 
may not be defined by elementary beam bending and twisting theory but 
yet whose effective restraint is amenable to evaluation. The deflec- 
tional and rotational stiffness provided by longitudinal stiffeners 
and full-depth webs is discussed and numerical exan^sles are given to 
illustrate the application of the charts to the design of wing structures . 


INTRODUCTION 


In current thin-wing construction, thick cover plates are often 
supported or stiffened by thinner gage members whose stiffness deter- 
mines the stability and strength of the cover plates . Experimental 
evidence is accumulating which indicates that, for certain longitudinal- 
stiffener proportions and for certain full -depth webs such as are found 
in multiweb construction, the supporting stiffness of these members must 
be carefully evaluated if reliable buckling stresses for the wing covers 
are to be calculated. The cross sections of these members often deform 
locally in a manner which permits deflection and rotation of the attached 
cover plate without developing the idealized beam or plate stiffness of 
the supporting member. For such construction, buckling stresses may 
not be predicted by stability criteria based upon beam stiffness quantities 
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such as El for tending and GJ for torsion, or by criteria based upon 
an assun5>tion of integral joints between cover plates and support members 
The desirability of relating plate stability to a stiffness parameter 
which defines the actual or effective stiffness provided by the support 
is therefore evident. 

The purpose of this paper is to present stability criteria for a 
number of supported-plate configurations frequently occurring in aircraft 
wing construction in which criteria the theoretical elastic buckling- 
stress coefficient is given as a function of the buckle length of the 
plate for a range of effective deflectional or torsional stiffnesses of 
the supports. In order to use the criteria for a particular supported- 
plate configuration, the effective restraint provided by the supports 
must be amenable to evalioation. A section of the paper is devoted to 
a discussion of procedures for evaluation of the deflectional and tor- 
sional stiffness provided by longitudinal stiffeners and full-depth webs, 
including the effects of cross-sectional distortion. Nxmerical examples 
are then given which illustrate this evaluation for practical design 
problems. The derivations of the stability criteria are Included in 
the appendixes . 


SYMBOLS 


b width of plate between intermediate supports 

A length of buckles 

B = A/b 

t thickness of plate 

x,y coordinate axes in length and width directions, 

respectively 

w deflection normal to plane of plate 

p nvimber of bays in width of plate 

q number of buckles occurring across width of plate 

n integer 

an Fourier coefficients 

N compressive load per unit width acting in x-direction 

(length direction) required to cause buckling 
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nondimens ional buckling-load coefficient, Nb^/n^D 
compressive stress 
critical compressive stress 
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Young's modulus of elasticity 
Poisson's ratio 

plate flexural stiffness per unit width. 
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deflectional stiffness per unit length of support, Ih/in.^ 

rotational stiffness of intermediate support (moment per 
unit length required to produce a rotation of one radian) 

rotational stiffness of edge support (moment per unit 
length required to produce a rotation of one radian) 


nondimensional deflectional restraint pareuneter 


ah 7 h 


nondimensional rotational restraint parameters 


6 


s,sl,sll,s^ 

c 

U^,U2,U5 


nondimensional rotational restraint parameter from 
reference 1 

plate edge rotational stiffnesses defined in reference 2 
plate carry-over factor defined in reference 2 
energies of deformation 


Vl 


work of applied stress 
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total potential energy of system 
energy parameter 
Lagrangian multipliers 

coefficients defining amplitude of support deflection 
cross-sectional area of stiffener 

moment of inertia of stiffener cross section about its 
own center of gravity 

modal coefficient affecting deflectional stiffness of 
longitudinal stiffener 

nondimensional bending stiffness parameter for stiffeners 
of sturdy cross section 

ratio of average stress in stiffener to average stress 
in plate 

Euler column load 
torsion constant 
shear modulus of elasticity 

torsion coefficient which takes into account bending 
stiffness 

polar moment of inertia 

amplitude of sinusoidally distributed lateral load 

lateral deflection of longitudinally compressed stiffener 
subjected to sinusoidal lateral load 

depth of web 

thickness of web 

Et^^ 

plate flexural stiffness per unit width of web, 

(1 - 1.2) 
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kyj buckling -stress coefficient of web 

z distance between center of gravity of stiffener and 

middle plane of plate 

p radius of gyration of stiffener about its centroid 

STATEMENT OF PROBLEM 


In figure 1 are shown portions of several wing cross sections in 
which the material carrying bending stress is mainly concentrated in 
the thick plates forming the wing contour. Running spanwise are a num- 
ber of lighter structural members in the form of longitudinal stiffeners 
and full-depth webs . In addition to carrying longitudinal stresses 
these members must resist cover-plate deflection and rotation at their 
respective locations by virtue of their stiffness. If the stiffness 
chaxacteristics of these members can be defined, the buckling stress 
for the construction can be calculated. 

In this analysis the assumption is made that longitudinal stiffeners 
and full-depth webs will provide a restraint to the attached cover plate 
which is proportional to the distortions of these support members. This 
condition is met if sinusoidally distributed normal loads or torsional 
moments on the supports are assumed to cause sinusoidally distributed 
distortions which are in phase with the loading. Thus support stiffness, 
which is the ratio of load intensity to distortion at any point, is a 
constant along the length of the support. With this support character- 
istic, the attached plate will buckle with deflections and rotations 
that are distributed sinusoidally in the length direction. 

A cross section of the cover-plate buckling modes considered most 
likely to occur are sketched at the right of each wing cross section in 
figure 1 ard are denoted cases 1 to 6. Cases 1, 2, and 3 primarily 
involve the deflectional stiffness characteristics of the support mem- 
bers, and cases 6 involve the torsional stiffness character- 

istics of the supports . For a given wing cross section, both modes of 
buckling should be investigated to determine which mode leads to the 
lower buckling stress . 

Cases 1 and 4 represent the buckling modes of a cover plate sup- 
ported by substantial shear webs with an intermediate spanwise member 
(shown as a longitudinal stiffener) centrally located between the webs. 
The shear webs are assumed to prevent deflection but may offer a tor- 
sional restraint to the cover plate. In case 1 the stability of the 
compressed plate was investigated for a range of deflectional stiffnesses 
of the intermediate support and in case k the torsional stiffness of the 
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supports was considered. Because the two lowest buckling modes are 
either symmetrical or antisymmetrical with respect to the spanwise center 
line of the plate, it is not necessary to consider both the deflectional 
and rotational stiffness of the support simultaneously. 

Cases 2 and 5 represent the most likely bucllLing modes for a cover 
plate with two equally spaced spanwise stiffening members of equal 
stiffness between shear webs. In case 2 the effect of support deflec- 
tional stiffness was investigated by assuming the torsional stiffness 
of the intermediate supports to be zero. The torsional stiffness of 
the intermediate supports was considered in case 5 with the assuD^jtion 
that the supports are capable of preventing plate deflection at their 
locations . 

Cases 3 and 6 represent the most likely buckling modes for a plate 
stabilized by many spanwise lines of support of identical stiffness . 

These supports may be full-depth webs, as indicated in figure 1, or 
longitudinal stiffeners. In case the deflectional stiffness of the 
supports was considered by assuming the support torsional stiffness to 
be zero. Torsional stiffness of the supports was considered in case 6 
in which the deflections along the supports are assumed to be zero. 

The loading and support conditions for the six cases considered 
are shown schematically in figure 2. The compression cover plate is 
represented by a uniformly compressed long flat plate which is simply 
supported at the loaded edges. The deflectional stiffness of the sup- 
ports is represented by an elastic spring whose stiffness per unit 
length is denoted by ijr. The stiffness may Include the flexibility 
of the tension cover of the wing in those constructions where this 
flexibility would have an effect on the stability of the compression 
cover. The parameter as defined in this paper is a generalization 

of the foundation modulus concept as used by Timoshenko for beams on 
an elastic foundation (ref. 3)* The support torsional stiffness param- 
eters are denoted by 7 and a. The parameter 7 is associated with 
the torsional stiffness of the nondeflecting shear webs ajid a is 
associated with the torsional stiffness of the intermediate supports . 
These two parameters are equivalent to the torsional stiffness param- 
eter i+SQ defined by Lundqulst and Stowell in reference 1. 

For each of the first three cases a stability criterion in closed 
form is derived by the Lagrangian multiplier method (ref. 4 ). For the 
last three cases a stability criterion is obtained by using the prin- 
ciples of moment distribution explained in reference 2. With these 
stability criteria, numerical calculations have been made and are pre- 
sented in design-chart form. 
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PRESENTATION OF STABILITY CRITERIA 


Cases 1, 2, and 3 »- The stability criteria for cases 1, 2, and 3 
which involve the deflectional stiffness of the intennediate supports 
are presented in appendix A as equations (AI 9 ), (A24), and (A28) . In 
these equations, the effective deflectional stiffness \[f of the sup- 
ports is contained in the nondimensional parameter and the 

effective torsional stiffness y provided along the shear webs is con- 
tained in the nondimensional parameter yb/rt^. Values of the param- 
eter xlfb^/jt^D may be detemiined from these equations as a function of 


the compressive buckling-stress coefficient k = and the ratio of 

buckle length to bay width A/b for assigned values of the torsional 


restraint parameter yb/n^D. 


Two sets of numerical calculations have been made by assigning the 
values 0 and “ to yb/n^Dj these values correspond to simple support 
and coiT 5 )lete fixity, respectively, along the shear webs. These numerical 
results are presented in tables I, II, and III. Cross plots of the 
values in the tables have been made to form design charts (figs. 3 to 7). 
From these charts, the combinations of \Jrb3/n4i)^ and A/b at which 
buckling is initiated, may be read. The cutoffs in figures 5 7 

define the values of \|fb5/jt^D at which general instability involving 
deflection of the cover and the supports changes to local buckling of 
the cover (no support deflection) in accordance with the assun 5 )tion 
made that the supports possess zero torsional stiffness . 

In order to use the charts for plates on particular types of sup- 
ports, the parameter \|rb5/jt^D for the support must be evaluated. For 
the ustoal type of support, such as a longitudinal stiffener, or a full- 

depth web, \|rb5/jt^D will be a function of the stresses in the support 
and the wave length of budding, as well as the physical characteristics 

of the support. A discussion of the evaluation of ilrb^/jt^D for longi- 
tudinal stiffeners and webs is given in the section entitled Effective 
Stiffness of Supports," and numerical examples illustrating the procedure 
are given in a subsequent section entitled "illustrative Examples." 

Cases 4 , 3 > and 6 .- For cases 4 , 6 , the cover is restrained 

by equally spaced nondeflecting supports of equal rotation stiffness a 
while the plate side edges are restrained by nondeflecting supports of 
equal rotational stiffness y. The stability criteria for these cases 
are given in appendix B as equations (B2), (b 6), and (BIO). Values of 
the rotational stiffness parameter ab/n^D required to develop a given 
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con5)ressive buckling-stress coefficient k = in the cover at a 

given ratio of buclie length to bay width 'h/lo may be determined from 
these equations for assigned values of the edge-restraint parameter /b/jt^D. 
As was done for the deflectional stiffness cases, numerical results are 
presented for /b/jt^ equal to 0 and The numerical results were 

obtained by using the stiffness tables of reference 5 have been 
plotted to form design charts (figs. 8 to 12). 

For a given design problem in which the supports have both deflec- 
tional and rotational stiffness, the buckling-stress coefficient obtained 
by considering the mode of buckling which involves the rotational stiff- 
ness of the supports must be compared with the coefficient obtained by 
considering the mode Involving primarily the deflectional stiffness of 
the supports. The lower of these two values defines the buckling stress 
for the configuration. The evaluation of the torsional stiffness of 
longitudinal stiffeners and full-depth webs is discussed in the next 
section. 


EFFECTIVE STIFFNESS OF SUPPORTS 


General design charts have been presented, which, with one reserva- 
tion, are independent of the medium providing restraint to the compres- 
sion plate. The reservation is that the supporting medium must be of 
such a type that sinusoidally distributed normal loads and torsional 
moments cause sinusoidally distributed distortions which are in phase 
with the loading. Such behavior is characteristic of beam stiffness, 
as provided by longitudinal stiffeners of sturdy cross section. The 
distortions of the unstiffened webs during buckling of a multiweb beam 
also appear to be distributed sinusoidally along the length of the beam, 
and the reactions of the attachment flange on the compression cover of 
the beam are assumed to be proportional to the distortions. 

A thorough analysis of the stiffness characteristics including 
effects of cross-sectional distortions and shear distortion of these 
two types of supports is beyond the scope of this paper; however, it 
is shown how these effects may be included in an evaluation of the 
stiffness parameters and ab/n^. 

Stiffness of longitudinal stiffeners .- The most common type of 
supporting medi\am for plates is the longitudinal stiffener which par- 
ticipates in carrying the compressive load. If the distortion charac- 
teristics of such a stiffener are defined by elementary beam bending 
theory, the deflection under a lateral load of amplitude g distributed 
sinusoidally over a length A is 
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where crAg is the end load carried hy the stiffener, and Igff 

moment of inertia of the stiffener cross section about an axis lying in 
a plane parallel to the attached plate. The stiffness of the stiffener, 
defined as the ratio of lateral load to deflection, then is 


= 



( 1 ) 


If the average stress 0 in the stiffener is proportional to the com- 
pressive buckling stress acting in the attached plate, ijr may be 
written as 


^ = 




ck 


bt 



or 


. (|)Y5^2" - ok ^ 4) 

VA/\bD bt ^2/ 


( 2 ) 


where c is the ratio of the average stress in the stiffener to the 
average stress in the plate. 

A recent theoretical analysis by Seide (ref. 6) shows that the 
effective moment of inertia of longitudinal stiffeners attached to one 
side of a uniformly compressed plate may be expressed as a correction 
to the moment of inertia of the stiffener about its own center of 
gravity Ig . In this form, equation (i<-) of reference 6 may be written 
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Eleff 

bD 



(3) 


In eqxiation (5) the quantity z is the distance between the center of 
gravity of the stiffener and the middle plane of the plate, and p is 
the radius of gyration of the stiffener. The modal coefficient Zp^ 

is a function of buckling mode and associated wave length. The varia- 
tion of Zpq with A/b taken from reference 6 (which is applicable 

when the plate side edges are simply supported and when Poisson's ratio 
is 1/3) is given in figure 15 . The subscript p denotes the number 
of bays in the width of the plate, euid q denotes the nimber of buckles 
across the width of the plate (q is equal to 1 for the cases considered 
in this paper). With Eleff /bD defined, equation (2) should give satis- 
factory values of the stiffness parameter i|rb^/rt^D for stiffeners of 
sturdy cross section; that is, stiffeners whose cross sectional and 
shearing distortions under load Introduce deflections which are small 
compared with the over-all deflection as a beam. 

In practical applications stiffeners are often formed from sheet, 
which necessitates a bend radius between the web of the stiffener and 
the flange for attachment to the sheet. For certain proportions, deflec- 
tion of the plate may be appreciably increased by the flexibility of 
the attachment flange between the rivet line and the web of the stiffener, 
and by shearing distortion in the stiffener. If the total deflection 6 
is assumed equal to 5p + 62 + 85 where 8p is the deflection due to 

bending of the stiffener as a beam, 82 is the deflection due to flexi- 
bility of the stiffener attachment flange, and 85 is the deflection 
due to shearing distortion in the stiffener, the effective stiffness 
may be written as 



2 ^ 


In nondimens ional form the effective stiffness is given by 
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where is given by the right-hand side of equation (2), 

must be evaluated either analytically or experimentally, and 

ilr^b^/n^D may be calculated. It is evident that if either \|t 2, or 

\|fj approaches zero, the effective stiffness of the stiffener approaches 

zero. Any other significant distortions can be included in a similar 
manner. 

The torsional restraint furnished a plate by a stiffener which 
undergoes no cross-sectional distortion when it twists is discussed in 
reference 7 * The expression for its stiffness (eq. (8) of ref. 7 
rewritten in the notation of the present paper) is 


a = s|(gJ + 2|eCbt - 



where the quantities J, Cbt, and Ip must be calculated with respect 
to an assumed axis of rotation. In nondimensional form, the stiffiless 
is 


^ = ftxVoj ^ ,2 b| 

n^D U/ \bD A^ b^D b^t/ 


( 5 ) 


Expressions similar to equation (5) should be derived for those stiffeners 
in which torsional moments applied to the stiffener attachment flange 
cause distortion of the cross section of the stiffener when it twists. 

Stiffness of full-depth webs .- Wlien the compression cover of a 
beam is supported by full-depth webs as in a multiweb beam, the effective 
stiffness of the webs in resisting sheet deflection and rotation at the 
sheet-web juncture must be evaluated. Reference 8, for exaii5)le, evaluates 
the effectiveness of integrally joined webs as torsional restraints on 
the cover of a multiweb beam. The assun^ition made in that einalysis is 
that the webs possess sufficient deflectional stiffness to form longi- 
tudinal nodes along the sheet-web juncture during buckling. The range 
of sheet and web proportions for which this assumption is valid, however, 
is not established. 

For built-up construction, the deflectional stiffness provided by 
an unstiffened web plate is influenced by the eccentricity of the con- 
nection between web and cover plates and by the state of stress existing 
in the webs of a beam under load. In particular, for channel-type webs 
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formed from sheet, appreciable distortions of the attachment flanges 
and lateral deflection of the web are produced by either depthwlse 
crushing or stretching forces. In accordance with the stiffness analysis 
for longitudinal stiffeners, the stiffness of the channel should be 
analyzed under the action of a depthwlse load applied sinusoidally along 
the length of the attachment flange in the presence of the stresses that 
exist in the web during beam bending. This procedure is illustrated by 
a numerical example in the next section. The outcome of such an analysis 
is influenced rather strongly by the assumed eccentricity of the applied 
load (with respect to the plane of the web) and by the degree of clamping 
that is assumed to be provided by the riveted connection between web 
attachment flanges and the cover plates . The importance of these fac- 
tors in calculating deflectional stiffnesses has been emphasized by 
Bijlaard and Johnston in a recent paper (ref. 9 )- 

With regard to the torsional restraint provided to the con5)ression 
cover by integrally joined webs, the restraint data presented in fig- 
ure 9 of reference 8 are analagous to equation (5) for the torsional 
stiffness of a stiffener; that is, the restraint coefficient e in 
figure 9 of reference 8 is a measure of the negative of the stiffness 
of a web subjected to a pure bending stress distribution as a function 
of buckle length. The relationship between the torsional stiffness 
parameter ab/jt^D of the present paper and the restraint coefficient e 
is 

gb e b % 

byj D 


When webs are not integrally joined to the cover, the stiffness of 
attachment should be taken into account when the parameter gb/jt^ 
calculated . 


( 6 ) 


the 

is 


ILLUSTEIATIVE EXAMPLES 


Some of the procedures outlined in the preceding section for calcu- 
lating the effective stiffness of supports will be illustrated in the 
solution of two conmon cover-plate stability problems. The first 
exaji5)le chosen considers the type of restraint offered by the webs of 
a multiweb structure and the second considers the effect of one sided 
longitudinal stiffeners on plate buckling. 

Buckling of a multiweb structure .- When the webs used in a multiweb 
wing are formed from sheet metal there is no assurance that the deflec- 
tional restraint provided to the beam covers by the formed channel webs 
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is sufficient to form longitudinal nodes along the web lines and thus 
to force buckling of the type denoted as case 6. The subsequent calcu- 
lations illustrate a simple procedure that may be used to investigate 
the possible occurrence of buckling in the mode denoted as case 5- The 
calculations apply to a multiweb beam tested in pure bending in the 
combined load testing machine of the Langley structures research 
laboratory. The beam had four identical channel webs (5 cells) and it 
is assumed that the analysis for a beam with an infinite number of cells 
can be applied. The physical dimensions of the beam are as follows: 


Cover width between webs, b, in 5*75 

Cover thickness, t, in 0.125 

Channel web depth, b^, in 2.08 

Channel web thickness, t^, in O.O 5 O 

Bend radius between web and attachment flange, in 0.20 

Diameter of web -attachment rivets, in 5/l8 

Pitch of web-attachment rivets, in 9/i8 

Distance between midplane of web and rivet line of 

attachment to cover, f (attachment flange assumed to be 
effectively clamped to cover at the inner edge of the rivet 

shanl^ when closely spaced rivets are used), in O .56 

Yoiing's modulus for the 75S-T6 aluminum alloy, psi .... 10. 5 X 10° 

Poisson's ratio for the material 0.553 


In accordance with the procedure outlined in the preceding section, 
the deflectional stiffness of the channel must be evaluated under the 
action of a sinusoidally distributed lateral load of amplitude g on 
the channel web in the presence of the existing bending stresses. This 
loading is shown in figure l4(a). The lateral loading is applied a 
distance f from the web plane, the distance at which the flange has 
been assumed to be completely fixed to the cover plate. In order to 
compute the deflection at a given cross section, the channel is idealized 
as in figure lk(b). The attachment flange is cut from the web and 
asBiuned to be flat and to be free of longitudinal conpression stress. 
(This stress is usually small in relation to the critical buckling 
stress of that portion of the flange between the rivet line and the 
web.) Also, since the buckle length is large conroared with the dis- 
tance f , the longitudinal bending stiffness of the flange will be neg- 
lected in computing the distortions at a given cross section. These 
distortions are shown in figure li^(c). The left-hand edge of the attach- 
ment flange is free but maintains a zero slope (to match the slope of 
the attached plate when buckling occurs in the mode denoted as case ^) , 
whereas the right-hand edge is supported against deflection and elasti- 
cally restrained against rotation by the torsional restraint a' . The 
restraint a' represents the resistance to rotation which the web 
offers the flange and is a function of both buckle length and the 
bending stress in the web. Because of the flexible corner radius that 
actually exists between the attachment flange and the web, the beam 
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cover is assumed to be eqiially free to deflect up or down with the 
attachment flange. Simple tension and compression loading tests on 
channels with comer radii verify this assumption. With these simpli- 
fying assumptions and boundaiy conditions, the deflection 6 at any 
cross section may be written 


&(x) 


£L 

12% 


1 + 4 


5l 

a'f 


1 + 


?5L 

a'f 


sin 5^ 

X 


The effective stiffness of the channel, defined as the ratio of lateral 
load to deflection, then is 


f^ , Drr 

1 + 4 ^ 

a'f 

or in nondlmensional form 



where e is the restraint parameter from figure 9 of reference 8 and 
is defined as 



NACA TN 2987 


15 


Substitution of the physical dimensions of the beam into equation ( 7 ) 
for /n^'D gives 


= 8.93 — -■ (8) 

rt^D 0.175e - k 


In order to obtain numerical values for \|rb5/jt^D, the quantity € must 
be read from figure 9 of reference 8. Values of e may be obtained 
which are compatible with the bending stress distribution in the beam 
if the stress in the extreme fiber of the web is assumed to be equal 
to the average stress in the beam covers and the lengths of the buckles 
in the webs and covers are equal. From these two conditions, the 
following equations result: 


— = = 1.80 ^ 


( 10 ) 


The lowest value of the buckling-stress coefficient k which simul- 
taneously satisfied equations (8), ( 9 ), and (lO) is the desired value 
and is found by a trial and error procedure. 

The first step in this procedure consists in determining by trial 
and error the value of k which satisfied equations (8), ( 9 ), and (lO) 
for an assumed value of A/b . Values of € are read from the curves 

of figure 9 of reference 8, and values of i|rb5/it^D are read from fig- 
ure 7 of this paper. This procedure is repeated for several assumed 
values of X/b. If this procediire is used, values of k equal to 3*35> 
3 . 25 , 3 * 26 , and 5*^7 sire found for assumed values of X/b equal to O. 7 , 
0.8, 0.9) and 1.0, respectively. The final step is to minimize k with 
respect to X/b . The minimum value of k for this mode of buckling 

(case 3 ) Is thus found to be 3.24 at ^ = O. 85 . 
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In order to determine the hue kling-s tress coefficient that would 
he obtained with buckling of the type denoted as case 6, figure 2 of 
reference 8 may he used to read the huckling-stress coefficient directly. 
The use of this direct-reading chart involves an assumption of an integral 
joint between the webs and the covers, and the indicated k value is 
4.1, which is considerably higher than the value 3*24 previously obtained. 

The actual experimental values of the buckling and failure stress 
for the exajiple beam were 


t^cr = 55>400 psi 
“failure = ^6,600 psi 


and the mode of buckling observed was that of the case 5* If* “the value 
k = 5*24 is substituted into the familiar buckling equation 


. /t\^ 

12(1 - p2)\b/ 


a buckling stress of 54,800 psi is obtained. 

Buckling of a plate with one-sided stiffeners .- In calculations of 
the buckling stress for plates with stiffeners attached to one side, the 
assumption is commonly made that the moment of inertia of the stiffeners 
may be calculated about the plane of attachment to the plate. The 
following example illustrates the procedure for obtaining the buckling 
stress of the plate-stiffener coiribination when this assxjmption is made 
and also the slight variation in the procedure which is entailed by 
using the expression from reference 6 for the effective moment of 
inertia of a one-sided stiffener. 

Consider the effect of two equally spaced longitudinal stiffeners 
of sturdy cross section on the stability of a long compressed plate which 
is simply supported along the unloaded edges and supported by deflec- 
tionally rigid transverse ribs at equal intervals along the length. 

Assume that the stiffeners and ribs offer no torsional restraints to 
the plate. The following physical dimensions are given: 
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Plate thickness, t, in O.188 

Plate width between stiffeners, b, in 4.70 

Rib spacing, in 50 

Cross sectional area of l/8-inch thick Z-stiffener, 

As, sq in 0.451 

Moment of inertia of stiffener about its centroid, 

0.205 

Radius of gyration of stiffener, p, in 0.686 

Moment of inertia of stiffener about plane of attachment 

to sheet, In.^ 0.524 

Distance between centroid of stiffener and centroid of 

plate, z, in O.956 

Young's modulus for the 75S-T6 aluminum alloy, psi . . . . IO.5 X 10^ 
Poisson's ratio for the material 0.555 


The deflectional stiffness of a longitudinal stiffener of sturdy 
cross section is given by (see eq. (2)) 

jj4jj \A/\bD bt ^ 2 / 


if the 

Eleff 

sheet. 


con^pressive stress in the plate and stiffener are equal. If 
is calculated about the plane of attachment of stiffener to 

is 

bD 


^^eff 

bD 


12[l - (0.55)^]](0.524)E 
(4.70)(0 .i88)5e 
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If the buckle length is taken to be the rib spacing, the numerical 
egression for \lrb5/n^D is 






0.1.31 (^' )\ 

4.70 X 0.188 v*T0/ 


= 0.1080 - 0.01195k 
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The value of k which satisfies this equation simultaneously with the 

A 30 

curves of figure 5 at — = ; is the desired value. By trial ajid 

t 1^.70 

error, a common solution is found at k = In order to verify 

that k = 5-55 Is the lowest huckling-stress coefficient, the analysis 
is repeated by assuming that two buckles occur between rib stations. 

In this particular example, this assun5)tion leads to a much higher 
value of k. 


The buckling-stress coefficient is now computed by assuming that 


^^eff 

bD 


is given by 


^^eff ^ 
bD bD 



1 + Z 

pq. bt 


In order that the modal coefficient Zpq may be read from the 

curves of figure 15, the buckle length must be assumed. The previous 
calculation indicated that the length of the buclcle is 50 inches and 
that it extends across the entire width of the plate. Thus, with p = 5> 

read from figure 15 is O.8O. 
is then 


Eleff ^ 12(1 - (0.533)i|(0.205)E 
4.70(0.i88)5e 



0.80 X 0.451 
U.70 X 0.188 


1 j A 50 
q = 1, and - = 

^ l<-.70 


the value of Zpq^ 
®Ieff 


From the data previously given, 


bD 


* (69.4)(2.594) = 166.0 
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With this value for the expression for \lfb 5 /jt^D is 

bD 


^ = 0.100 - 0.01195k 


By the use of figure 5 ^ "tbe value found for k is 5 - 25 . This value is 
about 8 percent lower than the value 3*55 obtained when the moment of 
inertia was rather arbitrarily chosen. (For other plate-stiffener 
combinations, the difference in the k values calculated by these two 
procedures can be either larger or smaller than the difference obtained 
in this numerical example.) 

The buckling stress for the stiffened plate is found by substituting 
the value of k = 5*25 into the buckling equation and a bucld.ing stress 
of 50)^00 psi for 75S-T6 aluminum alloy is predicted. 


CONCLUDING REMARKS 


Design charts have been presented which permit the evaluation of 
the compressive buckling stress of a long flat rectangular plate with 
various deflectional and rotational elastic line supports running 
lengthwise of the plate. In order to use the charts in a particular 
plate buckling problem, the restraint provided by supportipg elements 
such as angle and z-sections and full-depth webs like those used in 
multiweb wing construction must be evaluated. The evaluation of the 
stiffness of these supports has been discussed, and possible approaches 
for obtaining the required stiffnesses are presented. Numerical exan5)les 
have been Included to illustrate the type of procedures involved in 
computing buckling stresses . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 5 > 1955 - 
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APPENDIX A 

DERIVATION OF STABILITY CRITERIA FOR CASES 1 , 2, AND 5 


Although a set of stability criteria could be derived for the 
general case involving any number of lines of support either by solving 
the plate differential equation or by the Rayleigh-Ritz energy method, 
a desirable gain in simplicity is achieved by applying the energy method 
using Lagrangian multipliers (see ref. 4) to the Individual cases. The 
latter approach is shown in some detail for case 1, and variations in 
the method are indicated for cases 2 and 3* 

Case 1 .- An exact representation of the buckle pattern for case 1 
is given by the following series 


w = sin ^ y 

^ n=l,3,5 


an sin 


mty 

2b 


(Al) 


where the origin of the coordinate system lies along a side edge of the 
plate. The sinusoidal deflection along the plate center line may be 
written as 


w(x,b) = A sin 


nx 


(A2) 


and the slope along the side edges of the plate may be written as 


dw 


^(x,0) = B — sin — 
¥ 2b A 


(A3) 


Compatibility of equations (Al), (A2), and (A3) requires that 


n=l,5,5 ^ 



00 



n=l,3,5 


ann 


- B = 0 


ikk) 
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Using equation (A1) permits the so-called strain energy of tending 
stored in the buckled plate to be written as 


'A r'2fb 


Ui = 


'0 


^0 


dx dy = ^ jt^Xb 



+ n<= 




(A5) 

Using equation (A2) gives the energy stored in the deflectional restraint 
as 


U2 



(A6) 


and using equation (A5) gives the energy stored in the torsional 
restraints as 



(AT) 


The so-called external work done by the uniform compressive load N at 
buckling is 


Vl 



N 2b 
4 A 


00 



n=l,5,5 



(a8 ) 


The total potential energy may now be written as 


T = (U]^ + Ug + Uj - Vp) 
or 



an 


n=l,5,5 




\lfb5 






7b 




(A9) 
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where 




k = 




and 


\(fb^ 

T" 


and 


7b 

it^D 


are the nondimens ional 


deflectional stiffness and rotational stiffness parameters, respectively. 


The budding load is determined by the condition that the potential 
energy T' must be a minimum. Since the coefficients A and B depend 
upon the Fourier coefficients, a^j the expression to be minimized, is 


Q = T' 


-AW 







C oo 

Z 

= 1 , 5.5 



(AlO) 


where Ap and A2 are the Lagrangian multipliers . The potential 
energy T' is a minimum when 


= ^ = ^ = 

8an 8 a 8b ^2 


(All) 


8Q 

8an 




- Ap sin ^ - A211 = 0 


(n = 1, 3, 5, • • •“) 
(A12) 


_ p. 


+ Ap = 0 


(AI3) 


- P 7b 

SB" 


3^ + A 2 — 0 


(All^) 


8 q 

^p 


®-n ^ “ 

n=l,3,5 


A = 0 


(AI5) 
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c>i5»2 


00 


n=l, 5,5 


anil - B = 0 


(A 16 ) 


Eqiiatlons (A 12 ), (AI5), and (A 14 ) may be solved for A, «-nfl B, 

respectively, and these expressions substituted into the compatibility 
conditions (eqs. (AI5) and (A16)). Thia substitution results in the 
following two simultaneous homogeneous equations: 



(AI7) 



^1 

00 

n sin ^ 
2 



n2 


V = n 

1 

2 



2 

- k 

+ 

2 n=l, 5,5 , 

^ “2“ 

^1 n2 ' 

i + — B 


f ^ 

(A18) 


Each of the infinite sums in equations (AI7) and (A 18 ) are amenable to 
exact evaluation. Resolving the infinite series in equations (AI7) 
and (A18) into partial fractions yields the following forms: 


sin2 Ig. 


A. 


2 ^ 1 

n=l^,5 n2 - iL cp2 


00 

H 


n=l, 5,5 n2 + 02 


w 

Z 


n': 


8pcp^ 






z 


8 p02 


n=l, 5,5 n 2 - iL Cp2 n=l, 5,5 n 2 + iL. 02 


I 

I 
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00 



n=l,5,5 


n sin ^ 



k 


23 

/k 


00 (-1)^ sin ■ 

Z 

n=l, 2,5 n2 - 1 . <p 2 
-2 ^ 


00 

2 £ y 

v'k n-f7a,3 


sin as. 

n2 + .!L 92 
2 


where 


(p = - 1 

P 

0 = -y^i/k + 1 


By using equation ( 6 . 1 ^ 95 ) of reference 10 , the infinite series can he 
written in closed form. Thus, 



03 

Z - 

^ 3 ^ 5 / -1 


* -^(9 tan 9 + 9 tanh 0 ) 
P\/k 


o 


2 


4R 
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Substituting the closed forms of the infinite series into equations (AI 7 ) 
and (A18) and simplifying yield the following stability criterion: 


/ sin cp sinh 9\ S'/k 

jt2 y cos cp cosh 0 / 

«2d 


2/ 1 1 

«\cos cp cosh 0^ 


0 = 


-f— M 

rt\co8 cp cosh 0/ 


^ sin cp 
(p 

KCOB cp 


sinh 0\ 
cosh 0/ ^^5 


(AI 9 ) 


For given values of k, 3, and 7b/«%), the value of vlrb^/jt^D which 
causes the determinant to vanish is the desired value. When the side 
edges of the plate are simply supported, which is equivalent to 
-y't) 

setting = 0, the criterion reduces to 

rt%) 





sin cp sinh 9 

cp ^ 9 

cos cp cosh 0 


(A 20 ) 


In reference 11, a stability criterion is presented for the conqjres- 
sive buckling of simply supported plates with an arbitrary n\nriber of 
longitudinal stiffeners . When equation (A7 ) of reference 11 is applied 
to an infinitely long plate and written in the notation of the present 
paper it appears as 
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^ bD / 


4\/k 


which is equivalent to 


sin cp 
__ 


sinh 8 

0 


cos ^ - cos cp cos - cosh 0 
p p 


jAd 


4 ^ 


sin cp 

q> 


sinh 0 
0 


cos — - cos cp cos — - cosh 0 
P P 


(A21) 


when the stiffness of a stiffener is defined by elementary beam theory 
and the stresses in the plate and the stiffener are equal (see eq. (2)). 
Equation (A21) may be used for plates with simply supported side edges 
and with an arbitrary number of longitudinal supports . Equation (A20 ) 
may be obtained from equation (A21) by substitution of the proper 
values of p and q for case Ij that is, p = 2 and q * 1. 


For complete fixity of the side edges, 
criterion (eq. (A19)) reduces to 


= 00, the stability 

it^D 



sinh 0 
0 


icosh 0 


4\/k 


sin cp \ / 1 1 

cp I ^ v^cos q) ~ cosh 8/ 

cos cp/ q 2 sinh 8 ^2 sin qp 

9 ^ cp 

cosh 0 cos cp 


(A22) 


Solutions of these eqxaations and those to follow are facilitated by a 
tabulation of the functions cp, sin cp, cos cp, 0, sinh 0, cosh 0 
for appropriate values of the parameters k and p. These data sure 
provided in table IV. 
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Case 2 .- An exact representation of the deflection for case 2 is 
given by 


w 


00 

^ n=l,2,5 


an sin 


nrty 

5b 


(A25) 


Following the Batme procedure as for case 1^ two critera are obtained, 
one for symmetrical hucliLing and one for an anti symmetrical wave pattern. 
Calculations made by considering both modes of buckling indicated that, 
except for a very limited combination of values of k and A/b (k ^ ^ 
and A/b in the neighborhood of unity), buckling in a symmetrical mode 

requires the highest values of the stiffness parameter to 

achieve a given bucld_ing-stress coefficient k. Thus, for most practical 
problems, the criterion for symmetrical buckling only need be considered 
and is given in determinant form: 


/ Bin (p b 1]^_9 \ 

^[(p2 1 + g cos q> ^ q 2 1 + 2 cosh 8 Q _ + 

<2 1 1 + COB 9 i - cos (p 1 ^ 2 ' ^ 


sinh 9 


Bln <p 

JP y 


9 cosh e ®^'P|-oos<p/ 


= 0 (A2lv) 



/ 

Birth 9 

5 

f ® 

e 

X ^ 

iBlnh 6 

i - cosh 9 


sin qp 


COB 9 , 


/ 

sinh 9 

( 

8 

VI 

- cosh 9 



When the plate side edges are simply supported, the criterion reduces to 




sin cp sinh 9 

q) 9 

— - cos cp 1 - cosh 9 
2 2 


(A25) 


which is the same as equation (A21) for q = 1, P ~ 5- 
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For complete fixity of the side edges, the criterion is 


sin cp 
i - COS cp 


sinh 8 

i* - cosh 8 
2 




Binh 8 
— 0 


sin cp 
9 


^slnh ® ~ - cosh 0 b1ii ^ 2 " 


02 


1 -f 2 cosh 0 


slnh 9 
0 


1 + cosh e 1 _ goal, 0 


+ V 


2 1 -»• 2 cos cp 


sin cp 

qp 


1 + cos cp 1 . cos cp 


(as6) 


Case 3 »- For the plate with many lines of support running longi- 
tudinally (case 3)> the stability will not be influenced by the side- 
edge conditions. Correspondingly, the following function is used to 
describe the deflection surface 


w 


sinf 2 : 


n=0,2,4 


"n 


cos 


njty 

b 


(A2T) 


where the origin of coordinates is taken midway between any two lines of 
support. Physically the problem thus considered is the buckling of an 
infinitely wide plate column of length A restrained against deflection 
along continuous longitudinal lines which are equally spaced across the 
width of the plate. The stability criterion for this case is 




sin cp 
9 


sinh 9 
0 


1 - cos cp 1 - cosh 0 


(A28) 


which is the same as equation (A21) for 


q = 1, p = 00 . 
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APPENDIX B 

DERIVATIDN OF STABILITI CRITERIA FOR CASES k, AND 6 


A direct way of obtaining stability criteria for cases h, 6 

is by application of the principles of moment distribution to the 
stability of plates as described in reference 2. For a long plate 
supported along longitudinal lines by nondeflecting supports, the 
stability criterion is obtained by setting the sum of the stiffnesses 
of the members entering the joint at a given support equal to zero. 

The plate stiffnesses are denoted in reference 2 by the symbol S, 
with appropriate superscripts, eind the carry over factors are given 
by the symbol C , with appropriate superscripts . These symbols and 
their superscripts will be used as defined in reference 2. The support 
torsional stiffnesses a and 7 as defined in this paper have an 
absolute value four times as large as S . 

Case 4 .- For neutral stability, the sum of the plate stiffnesses 
and the support stiffness at the joint along the plate center line must 
equal zero. The sum of these stiffnesses is 

i a + 2 S^ = 0 (B 1 ) 

k 

If equation (l2) of reference 2 is used, equation (B1) may be written as 


1 

jj- a + 




1 - C2 



vhich can be put into the foiloving nondimens ional form: 


A 

i(^ D 

C^ 


Sllb 



n^D 


ab 


1 


(B2) 
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Solutions to equation (B2) may be readily obtained by using the tabulated 
values of S^^b/D and C given in reference 5. 

For the particular case of simple support along the plate side 
edges, 2!^ = 0, equation (B2) reduces to 

«2 d 


5^ 8 £Ib ,0 

,2i) i? ® 


(B5) 


With complete fixity of the side edges, 


7b 


00, equation (B2) reduces to 


_8_ S^^b 

ab jt2 D 


= 0 


or, making use of equation (I5) in reference 2 gives 


gb ^ 8 Sb _ 


(B4) 


With the aid of the tabulated values of S^^b/D and Sb/D given in 
reference 5j equations (B3) and (b 4) have been plotted as curves giving 


the bucld.ing-load coefficient 

constant values of ab/jt%D. 
and 9* 


k = as a function of A/b for 

These curves are presented as figures 


8 


Case 5 ‘- If the stiffnesses of the members meeting along one of 
the intermediate lines of support (fig. 2) is s\immed, the following 
equation for neutral stability is obtained 


i a + = 0 

4 


(B5) 
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With defined by equation (l2) of reference 2, equation (B5) inay be 

written as 


1 

k 


a + 



= 0 


which can be written in the nondlmensional form 


gb 


4 




1 - 


It 


2 D 


(b6) 


siib 

it2 7^ 
It 2 d 


This stability criterion is readily solved by using the tabulated values 
of S^^b/D, S^b/D, and C given in reference 5* 


When 


7 b 


is equal to zero, equation (b 6) reduces to 


gb 4 (s^h S^b'\ 

«2d «2V D D y 


= 0 


(bt) 


and when 


gb 
It 2 d 


00 




the stability criterion is 
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Equations (B7) and (b 8) have been plotted in figures 10 and 11 and are 

Jfb2 

presented as curves giving the huckling-load coefficient k = as 

functions of A/b for constant values of ab/rt^D. 

Case 6 .- For a plate with many longitudinal lines of support 
(case 6) f the condition that the stiffnesses at a joint must vanish 
for neutral stability is given by 


r a + 2SlV = 0 

k 


In nondimens ional form, equation (B9) nay be written as 


^ _8 S^b ^ Q 

jt2 D 


(BIO) 


With the aid of the tabulated values of S^^^b/D given in reference 5^ 
equation (BIO) has been plotted as curves giving the budding- load 

coefficient k = as a function of A/b for constant values of 

ab/rt^D. These curves are presented in figure 12. 
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TABLE I 

VALUES OF BEFLECTIONAL STIFFNESS PAEAMETER FOR CASE 1 





for k 












X/b 








0 

1 

2 

3 

4 

5 

6 

Sin^ily supported side edges; ^ 0 

0 A 

-19.89 

-17.85 

-15.67 

-13.32 

-10.71 

-7.690 

-5.928 

.5 

-10.19 

-8.526 

-6.680 

-4.534 

-1.827 

2.253 

12.26 

.6 

-5.899 

-4.497 

-2.861 

-.7810 

2.376 

9.977 


• T 

-3.721 

-2.517 

-1.050 

.9544 

4.558 

18.79 


.8 

-2.504 

-1.461 

-•1573 

1-709 

5.309 

25.10 


•9 

-1-773 

-.8685 

.2748 

1.927 

5.076 

18.52 


1.0 

-1.309 

-.5230 

.4676 

1.877 

4.411 

12.58 


1.2 

-.7887 

-.1916 

.5426 

1.527 

3.052 

6.177 

20.15 

1.6 

-.3826 

-.0197 

.3999 

.9039 

1.542 

2.411 

3.760 

1.8 

-.2948 

-.0034 

.3254 

.7062 

1.162 

1.733 

2.495 

2.0 

-.2384 

0 

.2658 

.5610 

.9052 

1.309 

1.807 

2*5 

-.1624 

-.0081 

.1571 

.3353 

.5293 

.7427 

.9803 

3.0 

-.1268 

-.0192 

.0937 

.2127 

.3387 

.4725 

.6157 

4.0 

-.0956 

-.0349 

.0275 

.0917 

.1580 

.2263 

.2970 

6 .0 

-.0754 

-.0488 

-.0215 

.0062 

.0342 

.0626 

.0^14 

8.0 

- .0688 

-.0542 

-.0589 

-.0235 

-.0080 

.0076 

.0233 

10.0 

-.0672 

- .0568 

-.0471 

-.0373 

-.0274 

- .0175 

-.0076 

Clamped side edges; “ 

t(2D 

0,4 

-19.895 

-17.85 

-15.67 

-13.33 

-10 .73 

-7.766 

-4.155 

.5 

-8.541 

-6.717 

-A -626 

-2.092 

1.338 

7.154 

.6 

-5.9169 

-4.5417 

-2.960 

-1.079 

1.417 

5.371 

14.86 

• 7 

-2.606 

-1.266 

.3898 

2.658 

6.424 

16.11 

.8 

-2.565 

-1.596 

-.4663 

.9286 

2.818 

5.828 

12.499 

• 9 

-1.04l 

-.0965 

1.0515 

2.555 

4.777 

8*865 

1.0 

-1.4087 

-.7233 

.0654 

1.004 

2.186 

3.810 

6.376 

1.2 

-.4208 

.1383 

.7765 

1.529 

2.455 

3.667 

1.4 


-.3039 

.1060 

.5594 

1.070 

1.659 

2.358 

1.6 


-.2555 

.0559 

.3924 

.7602 

1.168 

1.626 

2.0 

-.4133 

-.2262 

-.0304 

.175^ 

.5928 

.6254 

.8696 

2.5 


-.2251 

-.0997 

.0276 

.1593 

.2957 

.4374 

5.0 


-.2263 

-.l4l6 

-.0549 

.0338 

.1246 

.2177 

k ,0 

-.2796 

-.2330 

- .1859 

-.1381 

- .0897 

-.0408 

.0088 

6.0 

-.2604 

-.2398 

-.2190 

-.1982 

-.1771 

-.1560 

-.1347 

8.0 

-.2542 

-.2426 

-.2308 

-.2186 

-.2074 

-.1958 

-.1839 

10.0 

-.2506 

-.2438 

-.2365 

-.2289 

-.2215 

-.2159 

-.2065 






NACA TN 2987 


55 


TABLE II 

VALUES OF DEFLECTIONAL STIFFBESS PARAMETER FOR CASE 2 


A/b 



\lrb 5 /it 4 ) for k 











0 

1 

2 

3 

4 

5 

6 




Simply supported side edges 

W 

II 

0 



0.4 

-19.826 

-17.744 

-15.491 

-15.029 

-10.234 

-6.927 

-2.654 

.5 

-10 .050 







.6 

-5-707 

-4.1855 

-2.568 

-.0123 

3.511 

10.600 


.8 

-2.275 

-1.1166 

.3543 

2.545 

5.626 



1.0 

-1.094 

-.2287 

.8404 

2.265 

4.4295 



1.2 


•04579 

.8196 

1.790 

5.094 

19.46 


l.k 


.1265 

.6970 

1.377 

2.223 

4.82 


1.6 


.1417 

.5747 

1.0716 

1.656 



2.0 

-.1279 

.1218 

.5924 

.6882 

1.015 

1.381 

1.796 

2.5 


.0877 

.2567 

.4356 

.6261 

.8297 

1.048 

3.0 

- .04899 

.0624 

.1780 

.2982 

.4236 

.5545 

.6915 

4.0 


.0525 

.0964 

.1597 

.2284 

.2968 

.5667 

6.0 


. 00 ^ 

.0565 

.0648 

.0955 

.1221 

•1513 

8.0 


-.000436 

.0152 

.0310 

.0468 

.0628 

.0780 

10.0 

- .01465 

- .00466 

.00534 

.0154 

.0254 

.0556 

.0457 

Tb 

Clamped side edges; “ 

«2d 

0.4 

-19.623 

-17.745 

-15.495 

-15.035 

- 10.248 

-6.963 

-2.758 

•5 


- 8.548 

-6.527 

-5.980 

-.9679 

5.466 

12.90 

.6 

-5.5807 

-4.206 

-2 .418 

-.1385 

5.178 

10.118 


.8 

-2.1971 

-1.173 

.2171 

2.107 

5.371 



1.0 

- 1.0595 

-.3041 

.70^ 

2.0576 




1.2 

-.5948 

-.0534 

.70025 

1.650 

5.051 



1.4 

-.5772 

.0494 

.5919 

1.247 

2.122 



1.6 

-.2657 

.0686 

.4813 

.9603 

1-550 



2.0 

-.1602 

.0554 

.5145 

.5995 

.9209 

1.296 

1.767 

2.5 


.0268 

.1890 

.3614 

.5464 

•7475 

.9689 

5.0 


.00478 

.1159 

.2517 

.5550 

.4807 

.6158 

4.0 

-.07566 

-.0215 

.0399 

.1030 

.1675 

.2334 

.5015 

6.0 

- .06696 

-.0430 

-.0160 

.0115 

.0388 

.0666 

.0947 

8.0 


-.0515 

- .0358 

-.0207 

-.00540 

.00996 

.0245 

10.0 

-.06556 

-.0544 

-.0451 



-.0555 

- .0257 

1 

-.0159 

- .00619 
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TABLE IV 

VALUES OF FUNCTIOKS APPEARING IN THE STABILITY CRITERIA 


A/b 

9 

sin cp 

cos cp 

e 

slnh 9 

cosh 9 

k = 1 

0.4 

6 .08571 

219.331 

219.33 

9.2929 

5^30.57 

5^30.37 

.5 

4.44291 

42.5051 

42.517 

7.6953 

1099 .01 

1099.01 

.6 

3.31151 

15.6951 

13.731 

6.6251 

576.14 

576.14 

•7 

2.45821 

5-79911 

5-8846 

5.8516 

173.89 

175.90 

.8 

1.75621 

2.80881 

2.9815 

6.2686 

97.070 

97.075 

.9 

1.10581 

1.54201 

1.6756 

4.8115 

61.454 

61-462 

1.0 

0 

0 

1.0 

4.4429 

42.505 

42.517 

1.2 

1.1708 

.92106 

.58941 

3.8831 

24.277 

24.298 

1.4 

1.4192 

.98855 

.15102 

5.4764 

16.156 

16.187 

1.6 

1-5209 

.99875 

.04987 

5.1660 

11.835 

11.877 

2.0 

1.5708 

1.0 

0 

2.7207 

7.5626 

7.6277 

2.5 

1.5391 

.99950 

.05169 

2.3510 

5-2004 

5.2957 

5.0 

1.4810 

.99597 

.08967 

2-0944 

5.9987 

4.1218 

4.0 

1.3604 

.97795 

.20885 

1.7562 

2.8089 

2.9815 

6.0 

1.1708 

.92106 

.58941 

1.3853 

1.8729 

2.1231 

8.0 

1.0590 

.86190 

.50708 

1.1781 

1.4702 

1.7780 

10.0 

.9425 

.80905 

.58777 

1.0419 

1.2409 

1.5937 

k = 2 

0.4 

5.17601 

88.4851 

88.491 

9.8275 

9268.56 

9268.56 

•5 

5.40041 

14.9711 

15 .007 

8.2094 

1837.67 

1837.67 

.6 

2.05781 

3.77171 

5.9020 

7.1189 

617.54 

617.54 

•7 

.449921 

.465251 

1 . 1029 

6.3310 

280.86 

280.86 

.8 

1.4255 

.98914 

.14696 

5.7331 

154.46 

154.46 

.9 

1.8252 

.96851 

-.24973 

5.2619 

96.422 

96.427 

1.0 

2.0219 

.89996 

- .45595 

4.8813 

65.897 

65.905 

1.2 

2.1858 

.81677 

-.57696 

4.2995 

36.825 

56.858 

1.4 

2.2215 

.79577 

-.60558 

3.8737 

24.050 

24 .071 

1.6 

2.2064 

.80471 

-.59366 

5.5467 

17.335 

17.364 

2.0 

2.1240 

.850^ 

-.52541 

3.0735 

10.786 

10.832 

2.5 

2.0010 

.90888 

- .41705 

2.6762 

7.2305 

7.2993 

3.0 

1.8857 

.95082 

-.50972 

2.3977 

5 -iv 554 

5.5444 

4.0 

1.6949 

.99251 

-.12379 

2.0264 

5.7275 

3.8593 

6.0 

1.4525 

.99045 

.13786 

1.6126 

2.4082 

2.6076 

8.0 

1.2612 

.95246 

.30467 

1.5780 

1.8574 

2.1095 

10.0 

1.1589 

.90817 

.41859 

1.2225 

1.5506 

1.8451 
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■TABLE IV 

VALUES OF FUNCTIONS APPEARING IN THE STABILITY CRITERIA - Continued 


X/b 

9 

sin cp 

cos cp 

0 

sinh 9 

cosh 9 




k = 

3 



OA 

1^.55301 

58.8501 

58.862 

10.219 

15720 .4 

15720 .4 

•5 

2.29981 

4.93591 

5.0562 

8.5851 

2670.5 

2670.5 

.6 

1.0371 

.86095 

.50872 

7.4771 

885.56 

885 . 56 

•7 

2.0685 

.87868 

-.47741 

6.6756 

596.42 

596.42 

.8 

2.4587 

.64645 

-.76297 

6.0654 

215.35 

215.55 

.9 

2.6095 

.50735 

-.86174 

5.5838 

155.04 

153.04 

1.0 

2.6879 

.43829 

-.89885 

5.1927 

89.975 

89.981 

1.2 

2 . 7188 

.41051 

-.91195 

4.5954 

49.410 

49.421 

l.k 

2.6786 

.44665 

-.89472 

4.1528 

51.798 

51.814 

1.6 

2.6152 

.50415 

-.86562 

5.8151 

22.654 

22.657 

2.0 

2.4658 

.62551 

-.78021 

5.3188 

13.796 

15.852 

2.5 

2.2932 

.75022 

-.66119 

2.9012 

9.0705 

9.1255 

3-0 

2.1451 

.85957 

-.54325 

2.6067 

6.7402 

6.8140 

k.o 

1.9122 

.94229 

-.55481 

2.2115 

4.5099 

4.6195 

6.0 

1.6047 

.99945 

-.05390 

1.7675 

2.8421 

5.0129 

8.0 

1.4081 

.98679 

.16198 

1.5156 

2.1615 

2.5816 

10.0 

1.2692 

.95486 

.29666 

1 . 54^+7 

1.7882 

2.0488 


k 


0.4 

5 . 5124 i 

16.7501 

r — — 

16.779 

10.537 

18864.5 

18864 . 5 

.5 

0 

0 

1.0 

8.8858 

5614.54 

3614.54 

.6 

2.5416 

.71735 

-.69671 

7.77662 

1179.76 

1179.76 

•7 

2.8585 

.29847 

-.95442 

6.9528 

525.04 

525.04 

.8 

5 .0418 

.09963 

-.99502 

6.3521 

281.17 

281.17 

.9 

5.1221 

.01949 

-.99980 

5.8410 

172.06 

172.06 

1.0 

5.1416 

0 

-1.0 

5 . 44 i 4 

115.581 

115.585 

1.2 

3-0976 

.04598 

-.99902 

4.8273 

62.455 

62.441 

1.4 

5.0106 

.15062 

-.99145 

4.5755 

59.690 

59.705 

1.6 

2.9125 

.22729 

-.97383 

4.0240 

27.953 

27.971 

2.0 

2.7207 

.40857 

-.91272 

3.5124 

16.750 

16.779 

2.5 

2.5133 

.58776 

-.80903 

5.0781 

10.856 

10.882 

3.0 

2.5416 

.71735 

-.69671 

2.7706 

7.9528 

8.0154 

4.0 

2.0780 

.87411 

-.48573 

2.5562 

5.2280 

5.5228 

6.0 

1.7566 

.98629 

-.16505 

1.8879 

5.2270 

5.5784 

8.0 

1.5209 

.99876 

.049876 

1.6191 

2.4252 

2.6253 

10.0 

1.5694 

.97979 

.20004 

1.4597 

1.9912 

2.2282 
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TABLE IV 

VALUES OF FUNCTIONS APPEARING IN THE STABILITI CRITERIA - Concluded 


A/b 


Bin cp 

cos cp 

e 

slnh 6 

cosh 9 




k = 

5 



O.lf 

2.55191 

6.57681 

6.4547 

10.810 

24768.18 

24768.18 

•5 

2.1587 

.85210 

-.55462 

9.1442 

4680.05 

4680.05 

.6 

5.0604 

.08110 

-.99670 

8.0125 

1508.94 

1508.94 

•7 

5.5742 

-.25051 

-.97306 

7.1881 

661.798 

661.798 

.8 

5.4879 

-.55942 

-.94062 

6.5580 

352.45 

352.45 

.9 

5.5125 

-.56227 

-.95206 

6.0585 

215.87 

215.87 

1.0 

5.4928 

-.54405 

-.95895 

5.6514 

142.54 

142.55 

1.2 

5.5967 

-.25255 

-.96765 

5.0244 

76.057 

76.044 

l.k 

5.2754 

-.15541 

-•99105 

4.5606 

47.815 

47.826 

1.6 

5.1524 

- .01081 

-.999952 

4.2010 

53.569 

55.384 

2.0 

2.9270 

.21295 

-.97706 

3.6745 

19.702 

19.727 

2.5 

2.6925 

.45455 

-.90075 

5.2260 

12.570 

12.609 

5.0 

2.5019 

.59695 

-.80228 

2.9074 

9.1272 

9.1819 

4.0 

2.2157 

.80056 

-.59952 

2.4767 

5.9090 

5.9950 

6.0 

1.8450 

.96264 

-.27078 

1.9881 

3.5825 

3.7195 

8.0 

1.6158 

.99908 

-.042990 

1.7067 

2.6646 

2.8461 

10.0 

1.4520 

.99259 

.11852 

1.5184 

2.1729 

2.5920 




k = 

= 6 



0.4 

l.ll 64 i 

1.56521 

1.6906 

11.051 

51542.81 

51542.81 

•5 

2.9787 

.16217 

-.98676 

9.5718 

5872.07 

5872.07 

.6 

5.5885 

- .45217 

-.90178 

8.2285 

1875.12 

1875.12 

•7 

5.7940 

-.60709 

-.79461 

7.5945 

815.515 

815.515 

.8 

5.8468 

-.64818 

-.76147 

6.7558 

429.518 

429.518 

.9 

5.8511 

-.65615 

-.77155 

6.2487 

258.67 

258.67 

1.0 

5.7825 

-.59776 

-.80167 

5.8548 

171.00 

171 .00 

1.2 

5.6459 

- .48519 

-.87550 

5.1962 

90.291 

90.296 

1.4 

5.4975 

-.54844 

-.93752 

4.7227 

56.252 

56.241 

1.6 

5.5548 

-.21159 

-.97734 

4.5549 

38.925 

58.956 

2.0 

5.1017 

.059882 

-.99920 

5.8151 

22.680 

22.702 

2.5 

2.8445 

.29274 

-.95619 

5.5540 

14.291 

14.526 

5.0 

2.6585 

.48214 

-.87609 

5.0257 

10.280 

10.529 

4.0 

2.5296 

.72566 

-.68805 

2.5808 

6 .5660 

6.6417 

6 .0 

1.9578 

.95341 

-.55882 

2.0745 

5.9175 

4.0451 

8.0 

1.6954 

.99249 

-.12250 

1.7822 

2.8875 

5.0556 

10.0 

1.5228 

.99885 

.047978 

i 

1.5865 

2.5405 

2.5452 
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Figure 2.- Six cases for which sta 


Case 6 


criteria are presented 






NACA TN 2987 





^)X> 


Figure 3.- Stability curves for case 1 with simply supported side edges. 
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Figiore 4.- Stability cixrves for case 1 with clamped side edges. 
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Figure 8 .- Stability curves for case 4 with simply supported side edges. 
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Figure 9 *" Stability curves for case 4 with clamped side edges. 
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Figiare 10.- Stability curves for case 5 with simply supported side edges. 
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Figure 11 .- Stability curves for case 5 'wi'tb clamped side edges. 







Figure 12.- Stability curves for case 6. 
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Figure 13 •- Functions appearing in expression for effective flexural 
stiffness of stiffeners attached to one side of plate (from ref. 6) 




